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Oh' Abstract 






The pseudoperturbative shifted - / expansion technique PSLET 
[16-19] is generahzed for states with arbitrary number of nodal zeros. 
Bound- states energy eigenvalues for two truncated coulombic poten- 

•i-H . 

^ ■ tials are calculated using PSLET. In contrast with shifted large-N 

C^ • expansion technique, PSLET results compare excellently with those 

from direct numerical integration. 



1 Introduction 



Attractive truncated Coulomb potentials 

( h = 1,2,3 and a is a truncation parameter ) are of special physical 

interest. They avoid the singularity at r = ( the crux of divergence problems 
) and serve as models for many interesting physical phenomena [1-10]. For the 
case b = 1, Eq.(l) reads 

the eminent cutoff Coulomb potential. In quantum-field theory, it has been 
suggested that if gravitational interactions of elementary particles are taken 
into account, there would be a gravitational cut off of Coulomb interactions 
resulting in a finite theory of quantum field. Eq.(2) represents therefore 
a nonrelativistic version of this idea. It may, moreover, be considered as 
an approximation of the potential of a smeared charge rather than a point 
charge. When 6 = 2, Eq.(l) implies 

y(r) = -J^^. (3) 

often called the laser-dressed Coulomb potential. A model that has been 
found useful for the study of the spectrum of a laser-dressed hydrogen-like 
atoms when exposed to an intensive non-resonant laser-field [5-10]. It has 
been shown that under Kramers-Henneberger canonical transformations the 
potential of such atoms can be recast into (3) with the truncation parameter 
a is related to the strength of the irradiating laser field [7-9]. The effective 
potential for scattering by a uniform spherical charge distribution is well 



simulated by (3). Moreover, it is very similar to the modified potential of 
the nucleus of a muonic atom due to finiteness in its size [5,11-15]. 

As the Schrodinger equation for neither of the potentials is amenable 
to a general analytic solution, one has to retain perturbation techniques or 
numerical methods to analyze their bound states. Mehta and Patil [2] had 
investigated analytically the s-state energy for the potential (2); however, 
no numerical results were obtained . Intensive analysis had been carried 
out by De Meyer and Vanden Berghe [3] and Fernandeze [4] on the bound 
states of (2). The shifted 1/N expansion technique had been employed to 
calculate the energy eigenvalues of potential (3) [6]. David Singh et al [5] 
have employed a numerical method to calculate the energy eigenvalues for 
the potentials (2) and (3). As such, these potentials are good candidates 
to be analyzed through an analytical ( often semi-analytical ) technique to 
resolve their underlying physical aspects. 

Recently, we have introduced a pseudoperturbative shifted-/ ( / is the 
angular momentum quantum number ) expansion technique ( PSLET ) to 
solve for nodeless states of Schrodinger equation. It simply consists of using 
l/fas a pseudoperturbation parameter, where J = I — P, and /9 is a suitable 
shift. The shift f3 is vital for it removes the poles that would emerge, at 
lowest orbital states with / = 0, in our proposed expansion below. Our 
new analytical, often semi-analytical, methodical proposal PSLET has been 
successfully applied to quasi - relativistic harmonic oscillator [16], spiked 
harmonic oscillator [17], anharmonic oscillators [18], and two dimensional 
hydrogenic atom in an arbitrary magnetic field[19]. 

Encouraged by its satisfactory performance in handling nodeless states, 
we generalize PSLET recipe ( in section 2 ) for states with arbitrary number 
of nodal zeros, n^. > 0. In section 3 we apply PSLET to treat the potentials 
(2) and (3) and we compare the results obtained by PSLET with the exact 



ones. We conclude in section 4. 



2 Method 



The prescription of our technique starts with the radial part of the time- 
independent Schrodinger equation in h= m = 1 units, 

1 d^ l{l + 1) 



+ ^ — - + V(r) 



^„^,,(r) = E„^,i^„^,,(r). (4) 



where / is the angular momentum quantum number, n,. = 0, 1, ■ ■ ■ counts the 
nodal zeros and V{r) is an arbitrary spherically symmetric potential that 
supports bounded states. 

Most textbook perturbation techniques manipulate the potential term to 
introduce a perturbating expansion parameter. To the contrary, our method 
keeps the potential arbitrary to the condition of being well-behaved. We 
invest the centrifugal term to play this role. Thus, with I = I — j3 { f3 to he 
determined in the sequel ), Eq.(4) reads 

(5) 

Apparently, the natural limit of Eq.(5) is the large-/ limit. At that limit the 

centrifugal term dominates over the kinetic energy and the potential terms. 

This results in a semiclassical motion of the particle in an effective potential 

1 1 
Veff = — 2 ~'~7^^(^)' *here Q is a constant that scales the potential at large-/ 

limit and is set, for any specific choice of / and n^, equal to f' at the end of 
calculations[19, 20]. Hence, the motion is concentrated about the minimum of 
Veffi say tq. Consequently, a coordinate transformation through 

X = l^/'^{r -ro)/ro, (6) 



will be in point. It worths mentioning that the scaled coordinates, Eq.(6), 
has no effect on the energy eigenvalues, which are coordinates independent. 
It just facilitates the calculations of both the energy eigenvalues and eigen- 
functions. Performing the coordinate transformation, Eq.(6), Eq. (5) reads 



T" -^ 



i^ + f^'W)' 



^UrA^) = -fEnr,l^nrA^)^ 



(7) 



V{x{r)) 



T2 



I' + {2p + 1)1 + /3{/3 + 1) I 



X 



"■' V ^ 77- 



+ gV{x{r)) 



(8) 



Expansions about x = 0, i.e. r = tq, yeild 

1 ^_^(n + r 



X 



E(-i) 



■a;"r"/l 



^° [' ' vt, 



n=0 



(9) 



vixir)) = y: 



'rfn^\ i_roxrj_r,/ 



n=0 \ "'o 



n! 



(10) 



Apparently, the expansions in (9) and (10) center the problem at the point 
To and the derivatives, in effect, contain information not only at tq but also 
at any point on the axis, in accordance with Taylor's theorem. It is also 
convenient to expand E as 



n=-2 



(11) 



Equation (7) thus becomes 



' n -t't 



■^^ + f^W)' 



^nr,l{^) 



oo 



^ KA 



(«) /-(n+l) 



^nAx) 



(12) 



with 



^-fV{x{r)) = rll 



1 , V{ro) 
2rl Q 



Jl/2 



-X 



V {ro)rlx 

Q 



+ 



3 2 V^rpy^x"^ 
2^ + 2Q 
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n=l ^ 



+ -o^E 



n=3 



:-ii 
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(ir" 
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l~{n^2)/2 



+ ;5(;5+i)f^(_i)"(!i±ll:r-r("+2)/2+(M±i^ 

n=0 2 2 



(13) 



where the prime of V{ro) denotes derivative with respect to Tq. Equation 
(12) is exactly of the type of Schrodinger equation for one - dimensional an- 
harmonic oscillator 



■^^4"^^^+«" + ^w 



XnXx) — KirXrir 



X 



(14) 



where P{x) is a perturbation - like term and ^o is a constant. A simple com- 
parison between Eqs.(12), (13) and (14) implies 



^o = l 



1 rlVjr,) 

2 Q 



^ 2/3 + 1 ^ PiP 



2/ 



(15) 



+ 



1 rlVjro) 

2 Q 
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oo 
n=2 



(16) 



and 



\ = ^2 y- ^(n) J--{n+l) 



n=-2 



Equations (15) and (16) yield 



(17) 



E, 



(-2) 
rird 



1 ^V{ro 



2r2 Q 



(18) 



E, 



(-1) _ £ 

nr,l 2 
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(19) 



e(°) = - 

nr,l ^2 



2 '^ 



(20) 



4:!/ = A£Vr.^ ; n>l. 



(21) 



(-2) 



Here Tq is chosen to minimize E^^^ ^ , i. e 



rfE, 



(-2) 



72 7;.{-2) 



"'■'' and — -^ > 0, 
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(ir2 



(22) 



which in turn gives, with / = y/Q 



l-P = ^rlV'{T,). 



(23) 



Consequently, the second term in Eq.(13) vanishes and the first term adds 



a constant to the energy eigenvalues. It should be noted that energy term 
PeI^i is the energy of a particle moving under the effect of Ve/f. Hence, it 
is , roughly, the energy of a classical particle with angular momentum Lz=l 
executing circular motion of radius Tq in the potential V{ro). This term thus 
identifies the leading - order approximation, to all eigenvalues, as a classical 
approximation and the higher - order corrections as quantum fluctuations 
around the minimum r^, organized in inverse powers of /. 

The next leading correction to the energy series, lE~i, consists of a 
constant term and the exact eigenvalues of the unperturbed harmonic oscil- 
lator potential Q'^x'^/2. The shifting parameter f3 is determined by choos- 
ing IE~ I ={). This choice is physically motivated. It requires not only the 
agreements between PSLET eigenvalues and the exact known ones for the 
harmonic oscillator and Coulomb potentials but also between the eigenfunc- 
tions as well. Hence 



/9 



1 



[rir 



)VL 



(24) 



where 



VL 



qoV"{ro) 



Equation (13) thus becomes 



(25) 






/ 



^^V{x{r)) = ril 



1 , V{ro) 
2ry Q 



('"■)( ^\ I -'"-f^ 



^t;(")(x)/ 
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(26) 



where 



.(°)(a;) = if]V + ^, 



(27) 



t;W(a;) = ~{2f3 + l)x - 2x-^ + 



rlV"\r o) 
QQ 






(28) 



and for n > 2 



,W 



(x) = (-l)«(2/5+l)i!i±llx" + (-l)'^^^^^(n-iy"-2) 



+5„x"+2, 



(29) 



5„ = (_l)n(!i±^+ ^°"^'^ c^"+'V^(ro) 



2 g(n + 2)! cir^+2 



(30) 



Equation (12) thus becomes 



'2dx2 



+ J2 t;('^)r-"/2 



n=0 



^n.,/(a:) = r^ 



oo 



.n=l 



"^nrAx) 



(31) 



Up to this point, one would conclude that the above procedure is nothing 
but an animation of the eminent shifted large-N expansion ( SLNT ) [19,20]. 
However, because of the limited capabilities of SLNT in handling large -order 
corrections via the standard Rayleigh-Schrodinger perturbation theory, only 
low-order corrections have been reported, sacrificing in effect its preciseness. 
Therefore, one should seek for an alternative and proceed by setting the wave 
functions with any number of nodes as 



^n.,/(a:) = Fn^^i{x) exp{Un,,i{x)) 



(32) 



Eq.(31) is readily transformed into the following Riccati type equation 



Fl^,{x) + 2F;,,(x)f/„^/x) I + F„,,Kx){ - \[UI^A^) + (f/:^,,(x))2] 



(33) 
where primes denotes derivatives with respect to x.It is evident that (33) 
admits solutions of the form 



oo rir — l 



n=0 p=0 



(34) 



uUx) = E u!c\x)n^/' + E G(:)(x)r(«+^)/^ 



ra=0 



n=0 



(35) 



where 



f/(«) 



a; 



ra+l 



rir \"" ' A^ ^m,n,nr 

m=0 



X 



2m- 1 



) ^0,n,nr 



(36) 



n+l 



^nr V"^/ 2-^ ^m,n,r. 



J2m 



.X 



m=0 



(37) 



Clearly the nodal zeros of the wavefunctions are taken care of by Fn^^x). 
For nodeless states, 



F^A^) = 1 



(38) 



which reduces (33) to the problem described in our previous work for nodeless 
states[16-19]. To illustrate how our method works for nodal states, we will 



10 



treate the one-node state. Upon substituting equations (34)-(37) with n,. = 1 
into (33), it reads 

-KM {En=oUicK^)i-^^^' + E^=oG'(:)(x)r-("+i)/^) - iF;;,(x) = o 

(39) 
The solution of (39) then follows from the uniqueness of power series rep- 
resentation . Therefore equating the coefficients of same powers of / and x 
respectively, (of course, the other way around works equally well), one gets 

Di,o,i = -uJ, 4°; = 0, Ui'\x) = -ux, (40) 

<-^l,0,l — ! Oq,! — ) \^^) 

Co,o,i = -(2Ci,o,i + 2/3 + 1), (42) 

UJ 

D,,, = - (^ - B,] (43) 

to \ to J 

Bi Ai = ^ (|^2Ai + a,o,iCi,„,, - |(2/3 + 1), ) (44) 

Ei°l = 4 ("C' + 1) + "SCiAi - 5^ - ^) (45) 

.... and so on. Obviously, one can calculate the energy eigenvalue and eigen- 
functions from the knowledge of Cm,n,nriDm,n,nr and a^"^^ i'^ ^ hierarchical 
manner. Nevertheless, the procedure just described is suitable for a software 
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package such as MAPLE to determine the energy eigenvalue and eigenfunc- 
tion up to any order of the pseudoperturbation series. 

Although the energy series Eq.(ll), could appear divergent, or at best, 
asymptotic for small I, one can still calculate the eigen series to a very good 
accuracy by performing the sophisticated [N,M] Pade' approximation [22], 



p^im 



Po + PiA + --- + Pm/F 



IM 



(46) 



to the energy series, Eq(ll). The energy series, Eq(ll), is calculated up to 

Ei^li/l' by 



En.,l = P<? + i?l°i + ■ ■ ■ + i?£!,/P + 0(1/P), 

and with the P^{1/1) Pade' approximant it becomes 

E„,,[4,4] = /^<? + P,^(l//). 
Our technique is therefore well prescribed. 



(47) 



(48) 



3 Truncated Coulomb potentials 



In this section we consider the truncated Coulomb potentials, Eqs.(2) and 
(3). 

Substituting Eq.(2) into Eq.(25), one gets 



Q 



Itq + 3a; 



V ^0 + tt 
Eq.(23) together with Eq.(24) give then 

1 



/ + 



l + (2n^ + l)l 



Itq + 3a 
ro + a 



J/2 



{r-Q + a) 



(49) 



(50) 
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When substituting Eq.(3) into Eq.(25), one gets 



Eq.(23) in turn gives 



/ 



n = . 


rl + W 
rl + a^ 




l + (2n, + l) 


rl + 4^2 
rl + a^ _ 


'o 


(ro + a)3/4 



(51) 



(52) 



Eqs.(50) and (52) are explicit equations in r^. Clearly, closed - form solu- 
tions of Eqs.(49)-(51) for tq are hard to be found ( which is often the case 
). Thus, we use numerical methods to resolve the issue ( hence the notion 
that PSLET is often semianalytical ). Once ro is determined the coeffi- 
cients Cm,n,nr, Dm,n,nr, o.^p^l^ ^re determined in a sequential manner. Hence 
the eigenvalues, Eq.(47), and eigenfunctions Eqs.(34)-(37), are calculated in 
the same batch for each value of /, n^, and a. 

Tables 1 and 2 show PSLET results, Eq.(47), Epslet, for Is, 2s, 2p, 3p, 
3d, 4d and 4/ eigenstates of Eq.(l), covering week, intermadiate and strong 
ranges of the truncation parameter a. In addition, we display the Fade' 
approximants, Eq. (48), -E[4, 4], and the exact results, Eexact, obtained by 
direct numerical integration [5]. Comparing Epslet with Eexact, one notices 
the underlying relation between the accuracy of PSLET results and /, Ur and 
the truncation parameter, a. The accuracy of PSLET results increases with 
increasing / and/or n^. ( See figure (2) ). This is in accordance with our choice 
for the expansion parameter in Eq.(47), as l/l gets smaller as / and/or n^ 
increases. PSLET results show a good converging trend to the exact values 
as the truncation parameter gets larger ( See figure (1) ). In the strong range 
of a, PSLET results are almost exact. To resume the eigenenergy series, 
Eq. (47), Pade' approximant is calculated. Moreover, the stability noticed in 
different order Pade' approximants show that the results are accurate up to 
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ten digits at strong range of a, to the contrary of the of the results obtained 
in [5] where the accuracy is recorded up to eight digits in the same range. 

Table 3 displays PSLET results, Epslet, for some excited states of Eq.(3) 
along with the Pade' approximants, i?[4,4], of Epslet, the results obtained 
by SLNT[6], Eslnt and the exact results, Ef^^act, obtained by direct numer- 
ical integration [5] . Apparently, the accuracy of our results for this potential 
have similar behavior as in the previous case. When compared with the re- 
sults obtained by SLNT, our results show a better agreement with the exact 
ones. The difficulty in calculating high order corrections in SLNT through 
Rayleigh-Schrodinger perturbation theory results in loss in accuracy. PSLET 
makes it amenable to calculate high order corrections which improves the ac- 
curacy. 

Moreover, one notices that the results of PSLET are more accurate in the 
case of the truncated coulomb potential than that of the laser-dressed one. 
The reason behind this is that the truncated potential is more coulombic in 
nature which makes PSLET nearer to the exact results. 

4 Concluding remarks 

We have presented a generalization of our pseudoperturbative shifted - / 
expansion technique PSLET [16-19] to treat states with arbitrary number of 
nodal zeros, Ur > 0. Two truncated coulombic potentials have been treated 
via PSLET and very accurate energy eigenvalues are obtained. 

The outstanding features of the attendant PSLET are in order. 

It avoids troublesome questions as those pertaining to the nature of small 
parameter expansions, the trend of convergence to the exact numerical val- 
ues ( marked in tables 1-3 ), the utility in calculating the eigenvalues and 
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eigenfunctions in one batch to sufficiently higher orders, and the apphcabihty 
to a wide range of potentials. Moreover, beyond its promise as being quite 
handy, on computational and practical methodical sides, it offers a useful 
perturbation prescription where the zeroth- order approximation inherits a 
substantial amount of the total energy. 

Finally, the scope of PSLET applicability extends beyond the present 
truncated coulombic potentials. It could be applied to angular momentum 
states of multi - electron atoms [23-25] , relativistic and non relativistic quark 
- antiquark models [26], etc. We believe that the feature of our method in 
determining expressions for the bound-state wavefunctions makes it possible 
to describe electron transitions and multiphoton emission occurring in atomic 
systems in the presence of an intense laser field. It is therefore reasonable to 
re-examine such phenomena in the frame work of PSLET. 
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Tablel. Bound-state energies, in h= m = 1 units, of the potential 

V(r) = —- for the Is, 2s, 2p and 3p states. Epslet represents PSLET 

(r + a) 

results, Eq (47), -E44 is the [4,4] Fade' approximant, Eq. (48), and E^xact from 
direct numerical integration [5]. 





a 


Is 


2s 


2p 


3p 


—E PS LET 




0.387357746 


0.109481497 


0.117535370 


0.053309085 


—Eu 


0.1 


0.387922157 


0.109145059 


0.117535388 


0.053309210 


-^ exact 




0.38754365 


0.10950805 


0.11753535 


0.05330930 






0.180406651 


0.069577091 


0.082862488 


0.041787781 




1.0 


0.180368972 


0.069581801 


0.082862452 


0.041787655 






0.18036705 


0.06958066 


0.08286242 


0.04178766 






0.043439053 


0.024810349 


0.029446519 


0.018748152 




10 


0.043438645 


0.024810342 


0.029446516 


0.018748153 






0.04343872 


0.02481036 


0.02944652 


0.01879815 






0.012194732 


0.008579236 


0.009717589 


0.007237436 




50 


0.012194683 


0.008579235 


0.009717588 


0.007237436 






0.01219469 


0.00857924 


0.00971759 


0.00723744 






0.003653176 


0.002907080 


0.003169533 


0.002608781 




200 


0.003653168 


0.002907080 


0.003169533 


0.002608781 






0.00365317 


0.00290708 


0.00316953 


0.00260878 



Table2. Bound-state energies, in h= m = 1 units, of the potential 

V(r) = —- for the 3d, 4d and 4/ states. Epslet represents PSLET 

[r + a) 

results, Eq (47), i?44 is the [4,4] Fade' approximant, Eq. (48), and Eexact 
from direct numerical integration [5]. 





a 


M 


Ad 


4/ 


—E PS LET 

—Eu 

^ exact 


0.1 


0.054136568 
0.054136568 
0.05413657 


0.030648450 
0.030648449 
0.03064845 


0.030813599 
0.030813599 
0.03081360 




1.0 


0.045010006 
0.045010007 
0.04501001 


0.026625065 
0.026625059 
0.02662506 


0.027588160 
0.027588160 
0.02758816 




10 


0.021024302 
0.021024302 
0.02102430 


0.014373461 
0.014373461 
0.01437346 


0.015576600 
0.015576600 
0.01557660 




50 


0.007962796 
0.007962796 
0.00796280 


0.006153630 
0.006153630 
0.00615363 


0.006643882 
0.006643882 
0.00664388 




200 


0.002798562 
0.002798562 
0.00279856 


0.002353647 
0.002353647 
0.00235365 


0.002498272 
0.002498272 
0.00249827 
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Tables. Bound-state energies, in h= m = 1 units, of the potential 

V(r) = —-r-^ „-., ,„ for the "id, M and 4/ states. Epslet represents 

PSLET results, Eq (47), E44 is the [4,4] Fade' approximant, Eq. (48), Eslnt 
is from SLNT [6] and Eexact from direct numerical integration [5]. 





a 


2s 


3p 


Ad 


—E PS LET 




0.121234415 


0.055496250 


0.031244807 


— -E44 


0.1 


0.126937229 


0.055498046 


0.031244799 


—Eslnt 


— 


— 


— 


— 


^ exact 




0.12182090 


0.05549523 


0.03124480 






0.089679150 


0.052114869 


0.030781155 




1.0 


0.095048845 


0.052110018 


0.030781125 






0.089903 


0.052509 


0.030838 






0.09267933 


0.05206038 


0.03078150 






0.037161915 


0.028313951 


0.021412017 




10 


0.037154303 


0.028313599 


0.021412537 






0.037111 


0.028225 


0.021353 






0.03715440 


0.02831369 


0.02141257 






0.012450958 


0.010878999 


0.009476019 




50 


0.012450960 


0.010879000 


0.009476018 






0.016263 


0.010882 


0.009477 






0.01245096 


0.01087900 


0.00947602 






0.003923030 


0.003659489 


0.003410238 




200 


0.003923030 


0.003659489 
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(b) 



(c) 

Figure(l): The effect of the truncation parameter on the trend of conver- 
gence of Eq.(47) for Is-state of Eq.(2): (a)a = 0.1, (b) a = 50, (c) a = 200. 
Where the numbers on the horizontal axis represent the number of correc- 
tions added to the leading energy term of PSLET, Eq.(47), the vertical axis 
represents the energies ( in h= m = 1 units ) and the horizontal curve denotes 
the exact numerical results from [5]. 
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Figure(2): The effect of the angular momentum quantum number, /, and 
the radial quantum number, n^, on the convergence trend of Eq.(47) for (2) 
with a = 50: (a) Is-state , (b) 3p-state , (c) 4d-state. The figures prescription 
is similar to figure 1. 
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Figiire(l): The effect of the truncation parameter on the trend of convergence 
of Eq.(47) for Is-state of Eq.(2): {a)a = 0.1, (b) a = 50, (c) a = 200. Where 
the numbers on the horizontal axis represent the number of corrections added 
to the leading energy term of PSLET, Eq.(47), the vertical axis represents the 
energies ( in ft= m = 1 units ) and the horizontal curve denotes the exact 
numerical results fcom [5] . 
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